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ABSTRACT 

The possibility of measuring redshift space (RSD) distortions using photometric data 
have been recently highlighted. This effect complements and significantly alters the 
detectability of baryon acoustic oscillations (BAO) in photometric surveys. In this 
paper we present measurements of the angular correlation function of luminous red 
galaxies (LRGs) in the photometric catalog of the final data release (DR7) of the 
Sloan Digital Sky Survey II (SDSS). The sample compromise ^ 1.5 x 10^ LRGs dis- 
tributed in 0.45 < z < 0.65, with a characteristic photometric error of ~ 0.05. Our 
measured correlation centered at z — 0.55 is in very good agreement with predictions 
from standard ACDM in a broad range of angular scales, 0.5° < 6* < 6°. We find 
that the growth of structure can indeed be robustly measured, with errors matching 
expectations. The velocity growth rate is recovered as fas = 0.53 ± 0.42 when no 
prior is imposed on the growth factor and the background geometry follows a ACDM 
model with WMAP7+SNIa priors. This is compatible with the corresponding Gen- 
eral Relativity (GR) prediction fas = 0.45 for our fiducial cosmology. If we adopt a 
parameterization such that / — fl'J^{z), with 7 « 0.55 in GR, and combine our fas 
measurement with the corresponding ones from spectroscopic LRGs at lower redshifts 
we obtain 7 = 0.54 ± 0.17. In addition we find evidence for the presence of the baryon 
acoustic feature matching the amplitude, location and shape of ACDM predictions. 
The photometric BAO feature is detected with 98% confidence level at 2: = 0.55. 

Key words: data analysis - cosmological parameters - dark energy - large-scale 
structure of the universe 



1 INTRODUCTION 

The discovery of an accelerated cosmic expansion has be- 
come one of the biggest puzzles in modern cosmology over 
the last 10 years. Several scientific probes have been pro- 
posed to understand the nature of this acceleration. From 
"geometrical" tests based on measurements of the distance- 
redshift relation such as baryon acoustic oscillations (BAO) 
or Type la supernovae, to "growth" tests sensitive to the 
growth rate of perturbations such as redshift space distor- 
tions (RSD), weak lensing or cluster abundance. The success 
of these probes relies in the implementation of massive, and 
many times dedicated, observational campaigns that will 
scan a good fraction of the observable Universe. Some such 
surveys will base their science in galaxy redshifts derived 
spectroscopically, what provides accurate radial positions. 



* E-mail:martincrocce@gmail.com 



Others will instead measure redshift photometrically. This 
yields poorer determination of radial positions but allows to 
go deeper in redshift and have higher sampling rate. The 
later group involves the Dark Energy Survejj^ (DES), the 
Physics of the Accelerating Universe collaboratiorj^ (PAU) 
and the Panoramic Survey Telescope and Rapid Response 
Systeirj^ (PanStarrs) as well as proposals such as the Large 
Synoptic Survey Telescope]^ (LSST) and the imaging com- 
ponent of ES A/Euclid survey. 

Perhaps the most exciting results related to the large 
scale structure of the Universe to date have been obtained 
using spectroscopic data from surveys such as the two degree 
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field galaxy redshift survey (2dGRS, [CoUess et all ([200T|)) 



the Sloan Digital Sky Survey (SDSS, [York et al.| (|2000|)) or 
the recent WiggleZ Dark Energy survey fDrinkwater e t ah] 
1 2010 1). This is particularly so in regards to redshift space 
distortions and baryon acoustic oscillations, the two topics 
of this paper. 

Redshift space distortions arise because the receding 
speed of galaxies with respect to us is due not only to the 
Hubble expansion but also to their peculiar velocity. Hence, 
galaxy positions inferred with the Hubble law are modified 
with respect to their true positions depending on the lo- 
cal velocity field. At large scales the net effect results from 
the relative strength of the intrinsic clustering (the bias) 
and the amplitude of velocity flows set by the conserva- 
tion of mass through the growth rate of structure parameter 
/ — din D /din a (where D is the linear growth factor and 
a the cosmological scale factor). Photometric redshift errors 
are generally assumed to wash out these distortions. As re- 



cently discussed by Nock et al. (20101 and Crocce et al 
( 20Tl| ) this IS true to a good extent but does not remove 



the signal completely if one splits the data in redshift bins. 
Thus, RSD from photometric data could be a sensitive test 



for the growth of structure as a function of redshift (Ross 
|et al.|pMT| ). This can then be used to discriminate modi- 
fications of Einstein's gravity from dark energy models. 

In turn, BAO originate in the tight coupling between 
baryons and photons prior to recombination. At the time 
of decoupling their "last scattering" imprints a well defined 
comoving scale in the spatial distribution of baryons and 
matter of ~ 100 /i~^Mpc, characterized by a slight excess 
of pairs over random. This scale is today imprinted in the 
distribution of galaxies. Again, poor distance determination 
because of photometric redshift estimates wash out this ex- 
cess, at least in the radial direction. Photometric surveys 
should still be capable of detecting this signature in the an- 



gular distribution of galaxies (e.g. Blake & Bridle (2005l) 



In this paper we use the angular correlation function 
of luminous red galaxies (LRGs) in the imaging catalog of 
the final data release (DR7) of SDSS to address whether 
growth of structure can be robustly measured with photo- 
metric data despite several sources of systematic errors, low 
resolution photo-z and other unknowns. In parallel, we in- 
vestigate whether BAO can be observed in the clustering 
pattern of LRGs in concordance with model predictions af- 
fected by redshift space distortions. In a companion paper, 
[Carnero et al.| ( [2011[ ), we discuss the cosmological implica- 
tions of the presence of BAO in the clustering signal. 

These tests, that extend previous work using angular 

(l2007l) 



power spectrum (Padmanabhan et al. 



(20071; Thomas et al. (20111), may serve as a proof-of- 



Blake et al. 



concept for the potential of future, more accurate, photo- 
metric data to place interesting constraints into the nature 
of cosmic expansion, and/or provide valuable higher redshift 
leverage to complement spectroscopic measurements. 

This paper is organized as follows. In Sec. |2.1| we dis- 
cuss our data, including the selection of the galaxy sample, 
survey mask and photo-z. In Sec. [3] we describe our angu- 
lar correlation measurements including the error estimates 
and the impact of star contamination. Sec. |4] refers to our 
redshift space distortion analysis and the implications for 
the growth rate of structure. Sec. [5] is dedicated to discuss 



the evidence for the baryon acoustic feature in the measured 
correlation. Lastly, Sec. |6] contains our main conclusions. 



2 DATA 

2.1 Selection of the Galaxy Sample 

We perform a color based selection of LRGs in the photomet- 



ric catalog of the final SDSS II data release (DR7, Abazajian 
[etaL] ( [2059l )). We fol low two main steps. The first one, based 
on that published by Cabre et al. ( |2006[ ), aims at identifying 
the region in color-color space that is populated by high red- 
shift LRGs ( [Eisenstein et alT] ( |2001[ )) by selecting all those 
objects that verify 



(r-i) 



{g~r)> -0.72 {r-i) + 1.7, 



(1) 



where the variables g, r, i are model magnitudes corrected by 
extinction. The second step is a set of cuts that are intended 
to minimize the star contamination in the sample, 

17 < petror < 21, 

< (Jpetror < 0.5, 

< r - j < 2, 
< s - r < 3, 

22 < magso [mag/arcsec^] < 24.5, (2) 

where petror is petrosian magnitude (corrected by extinc- 
tion), apetror is the error on petror and magso is the sur- 
face brightness in magnitude petror at half-light radius rso 
(the radius containing 50% of petrosian flux), raagc,^ — 
petror + 2.5log{7vr^Q). These cuts yield a total of ~ 1.27x10'' 
objects with redshifts in the range ~ 0.4 to ~ 0.65. 

From the set of cuts in Eqs. ^ the first ones corre- 
sponds to our magnitude limits . The next ones ensure that 
colors correspond to a galaxy and effectively eliminates very 
few objects. Lastly the cut that is most effective, that in 
mag^o, with an upper bound on 24.5 to ensure well mea- 
sured galaxies and a lower bound on 22 to eliminate bright 
point-like objects (i.e. stars). We have found that eliminat- 
ing objects with mag^o > 24.5 mag/arcsec^ leads to no dif- 
ference in the clustering signal. However eliminating objects 
with mapso < 22 mag/arcsec^ reduces the amplitude of clus- 
tering at large scales by large factors even though they repre- 
sents a small percentage of the total sample. Hence, we next 
discuss the motivation for this cut in more detail (further 



evidence for this effect is given in Sec. 3.2 1 



The distribution of magso is given in Fig. IT] It is well 
concentrated around magso ~ 23 mag/arcsec but shows 
long tails due to objects contaminating the LRG sample. 
This contamination is more clearly depicted in the petror 
vs. mag^o diagram in Fig. |2] Top panel corresponds to 
our photometric sample and shows a different trend for 
magso < 22 mag/arcsec'^ and magso > 24.5 mag/arcsec'^, 
with the core of LRGs lying in between. Bottom panel shows 
the same diagram but for the SDSS DR7 spectroscopic sam- 
ple, after imposing the selection in Eq. ([T]). This panel nicely 
shows that all objects with mag^o < 22 mag/arcsec^ could 
be contaminated by stars. In addition. Fig. [3] shows a his- 
togram of number of objects per pixel (here the pixel size 
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Figure 1. Histogram of surface brightness (mag^o) for all the 
objects in our catalog 



is O.Oldeg^) as a function of galactic latitude and different 
redshift bins, including (solid) or excluding (dashed) galax- 
ies with low mag^o- Objects with low magso clearly concen- 
trate at low galactic latitudes introducing artificial density 
gradients towards the galactic plane (which then results in 
large density fluctuations at large scales). There is a slight 
gradient residual after imposing the cut in magso, which we 
avoid but leaving out galactic latitudes (denoted b) lower 
than 25°. We require, 

b>25°. (3) 
This yields a reduction of ~ 3% of the SDSS area used. 



2.2 Angular Mask 

We built the angul ar mask using a Healpix pixelization 
I Gorski et al] ( |2005[ ))f]over the entire sky with Nside = 512 
that yields a pixel size of ~ 0.01 deg^. 

We then eliminate from the analysis those pixels where 
the geometric acceptance of the survey is compatible with 
bad or no measurement by imposing a minimum number of 
galaxies per pixel {N /pixel) of 15. We notice that in order 
to build the mask we use all the objects in the photometric 
catalog (i.e. not limiting by petror < 21) because the den- 
sity of LRGs is very low to allow a robust and well pixelized 
mask construction. In addition we look at galactic extinction 
and magnitude errors maps in order to mask badly observed 
regions. Figure [4] shows the distribution of errors in r-band 
magnitude averaged in every pixel (i.e. mean error per pixel) 
in the top panel, and in galactic extinction in the bottom 
panel. There is a clear correlation between these two quan- 
tities in regions of high extinction. Hence we suppress from 
our mask pixels with bad mean error rather than applying 
the cut directly to the LRG selection, as this would imply 
introducing artificially low density regions and correspond- 
ing systematic effects. In summary we discard pixels with 
extinction higher than 0.2 and mean error higher than 0.3. 
We have checked that using different pixelization sizes for 
the mask, as well as diflerent levels of acceptance of a pixel 
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Figure 2. Surface brightness (magso) vs. Petrosian r appar- 
ent magnitude. Top panel corresponds to objects in our pho- 
tometric LRG catalog, bottom panel to objects classified as 
stars in the SDSS DR7 spectroscopic sample that verifies the 
same selection as the top panel. Top panel shows some dis- 
tinctive trends for magso < 22 mag/arcsec^ and for mag^o > 
24.5 mag/arcsec^, which contaminate our sample and can modify 
the clustering signal. Bottom panels makes clear that the region 
of magso < 22 mag/arcsec^ is populated by stars. The region with 
^agso > 24.5 mag/arcsec-^ correspond to badly measured galax- 
ies (not LRGs) and have no impact in our clustering analysis. 



into the mask (varying the threshold in extinction, mean er- 
ror and A''/pixel) does not change appreciably the measured 
angular correlation. 

The resulting angular mask is depicted in left panel of 
Fig. |5]in spherical equatorial coordinates, with right ascen- 
sion (ra) and declination (dec) along the x — y axis respec- 
tively. It spans from ~ 110° to 260° in ra and 75° of dec 
almost fully in the norther hemisphere. The vertical band 
at ra ~ 172° is due to the photo-z used in this analysis 
(C. Cunha, private communication), that is described in the 
next section. Notice that we only considered the largest con- 
tiguous area of the survey, discarding stripes 76, 82 and 86 
that contribute only a small fraction of the total SDSS area. 

This mask covers a total area of 7136 square degrees. 
We have verified that this angular mask is valid in the full 
redshift range used in this paper, and is therefore used in all 
the analysis presented here. 

Before moving onto the next section we include in the 
right panel of Fig. [5] the resulting galaxy density map in the 
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Figure 3. Histogram of number of galaxies per pixel (pixel size of 
0.01 sq deg) as a function of galactic latitude for different slices 
in redshift as indicated in the figure. In solid line we plot the 
histograms when we select galaxies with 18 < magso < 24.5, 
while in dashed-dot we cut 22 < magso < 24.5. Galaxies with 
low magso are contaminating low galactic latitudes {b < 25°). 




same coordinates as the mask and with the same healpix 
pixelization {Ngrid = 512). As expected, it looks homoge- 
neous over the whole area. 



2.3 Photo-z and Redshifts Distributions 

In this paper we us e the value added photometric catalog of 
( |2009| nthat is based on the Photoz2 tables o f 



Cunha et al. 

the full SDSS DR7 sample presented in |Oyaizu et~al] ( |2008[ ) 
This value added catalog was built using and extending 
the weighting method technique of Lima et al. (20081. As 



discussed in |Lima et al.| ( |2008| ) and |Cunha et al.||2Ci09| ) the 
technique aims at estimating the redshift distribution for a 
photometric galaxy sample (or selected subsamples) rather 
than estimating individual galaxy redshifts. Hence, as an 
added value the catalog provides accurate estimates of the 
redshift probability distribution, p{z), of each galaxy. 

We applied photo-z quality cuts to the catalog in or- 
der to remove badly defined p{z) (e.g. double or multi- 
ple peaked distributions that can represent outliers) as 
well as very broad ones (that can be interpreted as galax- 
ies with bad photo-z). To this end we impose two cuts, 
\zpeak — Zmean] < 0.05 and CTz < 0.1, where Zpeak is the 
peak of the distribution, Zmean is computed as J zp{z)dz 
and (Tz = J{z ~ Zmean)^ p{z) dz^ The first cut eliminates 
roughly ~ 11% of objects and the second ~ 9%. Imposed 
together these cuts reduce the sample by ~ 16%. These 
threshold values in and \zpeak — Zmean\ were obtained 
by identifying the tails in the distribution of values for these 
quantites in the full catalog. One can of course be more con- 
servative and impose more stringent cuts but at the expense 
of biasing the population towards brighter objects (that tipi- 
cally have better photo-z) or introduce shot-noise error due 



^ Available at http: / /www. sdss.org/dr7/products/value_added' 
* The catalog provides p{z) in 100 bins between z = 0.03 and 
1.47, hence these integrals are sums over these 100 bins. 




Figure 4. Top Panel shows a map of mean error in r-band mag- 
nitude per (Healpix) pixel. Bottom Panel is the distribution of 
galactic extinction. Noticeably the mean error is larger in zones 
of higher extinction (see text for details). 



to large decrements in the number of galaxies per bin. Our 
results, in terms of the of the best-fit models that we 
obtain, are robust and stable in front of the photo-z quality 
cut. 

In this paper we split the galaxy sample into redshift 
bins according to whether the maximum of p[z) lies in the 
bin or not. That is, we identify the maximum of p{z) as 
the photometric estimate of the true redshift [zphot] and 
do top-hat bins in photometric redshift. In turn, one of the 
most important ingredients in order to interpret the galaxy 
clustering signal is a robust estimate of the distribution in 
true (spectroscopic) redshift, N{z), of all the galaxies in each 
bin. Cunha et al. ( 2009 1 discuss several methods to obtain 



N{z) and shows, using both mock SDSS catalogs and train- 
ning spectroscopic subsamples, that their best and almost 
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Figure 5. Left Panel: The mask used in tliis analysis depicted with the Healpix mollview projection routine. The white regions are 
excluded from the analysis. Different grey levels display the 81 JK zones used as one of two methods to estimate the errors. The vertical 
band is due to the photo-z used in this analysis. Right Panel: Angular density map of the galaxy distribution in the photometric bin 
[0.5 - 0.6]. 



unbiased estimate is provided by the weighted sum of the 



training distribution (see also Lima et al.| ( ,2008J ) , which is 
equivalent to sum the p{z) distributions of all the galaxies 
in the photometric bin, 



(4) 



To test the accuracy in this determination of N{z) we 
have selected all those galaxies in our sample that are also 
included in the 2SLAQ spectroscopic catalog (~ 6000 ob- 
jects) and computed their distribution of true redshifts as 
well as N(z) according to Eq. (j4|. These two distributions 
are remarkably similar as shown in Fig. [6] A Gaussian fit 
to each of them shows that their peak differs by less than 
1% and their width by less than 9% 

This difference is in 

perfect agreement with the intrinsic scatter in true redshift 
distributions obtained from different photo-z codes (e.g. see 



Table Al in Thomas et al. (2011 1). Notice that we can not 
use 2SLAQ to estimate N{z) for our complete catalog since 
our LRG selection is different from that in 2SLAQ (in par- 
ticular the magnitude cuts). Nonetheless the previous study 
shows the degree of unknown in the red shift distribution. 
Hence we will use Eq. Q to estimate N{z) for our red shift 
bins and will discuss in Sec. |4.3| how our results vary when 
the width and/or peak change by 9% and 1% respectively. 

In this way, our selected sample of LRGs have a dis- 
tribution in true red shifts that peaks at z ~ 0.5 and ex- 
tends roughly from ~ 0.4 to ~ 0.65. For our analysis we will 



^ 2dF-SDSS LRG and Quasar survey, a stripe close to 0' 
nation within the imagini 
10 Peak and 



decli- 

Cannon et al.|||2006|) 



area of DR7 

idth are defined as the mean and standard devia- 
tion of the best-fit Gaussian distribution. 



From PDF addition 
Direct reconstruction 




Figure 6. Direct reconstruction of the true redshift distribu- 
tion for the spectroscopic sub-sample of our photometric cata- 
log against the estimate for the same distribution using individ- 
ual galaxy redshift probability distribution functions (PDF). The 
spectroscopic sub-sample was constructed from the 2SLAQ cata- 
log, Cannon et al. (2006), while PDF's are provided as a part of 
the DR7 SDSS value added catalog of Cunha et al. (2009). 



mostly refer to a single top-hat photometric redshifts bin in 
the range [0.5 — 0.6]. Figure [t] shows the true distribution of 
galaxies in this bin. The number of objects in this bin, after 
the photo-z quality cut, is 664870. Notice that the bin width 
is slightly larger than our typical photometric error at this 
redshift (a^ ~ 0.05) hence choosing a narrower bin would 
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yield almost the same distribution of galaxies but at the ex- 
pense of increasing the noise in the measurement of w{0) 
due to smaller number of particles in the bin (see below). 

We also consider three narrower redshift bins, with 
photoz ranges [0.45 - 0.5], [0.5 - 0.55] and [0.55 - 0.6]. The 
redshift distribution for these cases are shown in Fig. m They 
are clearly highly correlated and not narrower than the wider 
top-hat bin discussed above. The number of LRGs in these 
3 bins are 451753, 317882 and 346988 respectively. In addi- 
tion to an increase of shot-noise and overlap, one expects the 
estimation of N{z) to be not so robust for a bin narrower 
than the intrinsic photo-z (a possible evidence for this is dis- 
cussed in Sec. |4.3| |. These are the reasons why we decided to 
concentrate in a single redshift bin, and repeat our analysis 
in these 3 bins as consistency checks. 

Redshifts bins lower than z — 0.45 and higher than 
= 0.6 do not have enough number of LRGs to obtain pre- 
cise measurements of the angular correlation (we find 52845 
LRGs in the photometric range [0.4 — 0.45] and 30412 in 
[0.6 — 0.65]). At those extreme bins our measurements also 
become too sensitive to ours cuts (e.g. in galactic latitude 
and/or photo-z quality). This might be due to various rea- 
sons, for example, large magnitude errors that correlate with 
galactic latitude and lead to large photo-z errors. 

But perhaps the most worrisome issue happens at z > 
0.6 where we find a large extra-power over a broad range of 
large angular scales (already for 6 > 1°). Excess of poweip^ 
on large scales has already been found and discussed in dif- 
ferent LRG selections based on SDSS DR5 in Sawangwiti 
et al.| ( |2009|) and SDSS DR 7 i n [Thomas et al. | (2011 , 2oTo| 
(see aIso |Blake et al.| ( |2007| and|Padmanabhan ct al. (2007)). 



Nonetheless, we only encounter this problem in the redshift 
bin [0.6 — 0.65]. We defer a discussion of possible reasons 
for Appendix and proceed to discard these bins from our 
study hereafter. 



3 CORRELATION FUNCTIONS 

The angular correlation function measurements were per- 
for med s tarting from Healpix angular maps as described in 

{Nside 



Sec. |2.2| {Naide = 512, pixel size of ~ 0.01 deg ) and using 
a standard pixel estimator (Barriga & Gaztanaga (20021, 
[Eriksen et al] ( |2004l )) 

1 



(5) 



where Sq = Ng^i/Ngai — 1 is the fluctuation in number of 
galaxies in the i-th pixel with respect to the mean in the 
angular map, pixels i and j are separated by an angle 6 
and Npairs{6) is the corresponding number of pixel pairs. 
Pixels were weighted by or 1 according to the angular mask 
discussed in Sec. |2.2[ We have also implemented a standard 



Landy & Szalay ( Landy fc Szalayl|1993 1 estimator, and the 
resulting measured correlations were within 1% of that from 
Eq. ^. 

The measured correlations in the three bins of width 



We note that the exact meaning of excess is only loosely de- 
fined in the literature, in general is taken a roughly more than 
two sigma difference between model and measurements 
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Figure 7. True (spectroscopic) redshift distribution for the bin 
0.5 — 0.6 resulting from sum of the individual redshift probability 
distributions. A fit to a Gaussian function (shown in solid black) 
yields a media oi fi = 0.549 and standard deviation a = 0.062. 
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Figure 8. True redshift distribution for a set of "narrow" bins of 
width similar to the typical photometric error (Az = 0.05). 



0.05 as well as in the bin 0.1 are shown in Fig. [9] Error 
bars displayed in this figure were obtained using jack-knife 
resampling. In what follows we discuss our different error 
estimates. 



3.1 Error Estimates 

We estimate the error and covariance between angular bins 
in the measured correlation function using two independent 
methods. 

One method was implementing the standard jack-knife 
resampling technique, that to date has been widely used in 
clustering analysis with correlation functions (for a summary 
see |Cabri et al.| ( |2007[ ) and |Norberg et al.| ( [20091 ) ^nd refer- 
ences therein). To this end we divided the angular mask in 
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Figure 9. The measured correlation function in the different top- 
hat bins in photometric redshift. In this paper we focus on the 
bin [0.5-0.6] but test the robustness of our results throughout all 
the redshift range. 



81 jack-knife (JK) zones of similar area (~ 90deg^ each) and 
shape. These zones are shown in grey levels in Fig. [5] One 
then takes different realizations to be all the sampled area 
except from one JK zone at a time. The covariance is then 
computed from the dispersion among the measurements of 
w{8) in the Njk ~ 81 resulting realizations, 



Gov J 



Nj 



Nj 



Njk 



w{e)) 



y.{w%{e')-w{e')) 



where w{9) corresponds to the full area, 



„{») 



(6) 

to the i-th 



realization and the factor Njk — 1 corrects from the fact 
that realizations are not independent. The positive aspect of 
the jack-knife resampling is that this estimate is build out of 
the data itself and hence it can account for systematic effects 
difficult to capture otherwise. On the negative side is the 
fact that different realizations might share a large fraction 
of area. In principle this is accounted for by the rescaling of 
the covariance in Eq. ([6|, but this JK resampling may not 
reflect the true underlying statistical variance, particularly 
if the JK zones are too small or irregular 

It is then desirable not to rely only on the jack-knife 
estimator. Therefore we have also calculated the error and 
covariance matrix using the analytical approach discussed 
in Crocce et al. ( 2011[ ) in which. 



9 9/ -I- 1 

7^E(i^A(.)p.(.')(c. + iM)^ 



(7) 



where x = cos(S) and Ci is the analytical angular power 
spectrum for the redshift distribution of interest. This es- 
timate have been extensively compared against ensemble 
errors drawn from simulated photometric surveys assum- 
ing different binnings of the data, survey depth, underlying 
photo-z, shot-noise contribution and more, see [Crocce et al.| 
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Figure 10. Relative error in the measured correlation function. 
Symbols correspond to Jack-knife estimates using 81 zones while 
solid line to the theoretical expressions in Crocce et al. (2011). In 
the later case the ratio is against the best-fit to(S) model. 



( ]2011 1. In turn, Ross et al. (20111 recently showed that the 



distribution of best-fit values for the bias and growth rate of 
structure recovered in survey mocks agreed very well with 
the errors obtained when using Eq. Q. Furthermore, these 
studies used an angular mask comparable to the simple ge- 
ometry treated in this paper. 

To compute the angular spectra in Eq. ^ we assume 
a ACDM cosmology with WMAP7 parameters ( |Komats"u 
|etaI] ( |20TT| )) and use the redshift distributions in Figs.[7j|8| 
In turn, initial values for large-scale bias b and growth rate / 
are obtained from a minimization using jack-knife errors 
(see Sec|4|. Provided with the full Ci spectra we then com- 
pute the covariance matrix in Eq. ([7| with fsky = 0.1682. 

Figure 10 shows the relative error Aw/w in the mea- 
surement of the angular correlation for our main case bin 
in photo-z range [0.5-0.6] (where Aw = Cov(6l, S)^/^). 
The agreement with the relative error recovered from 
the jack-knife technique is remarkable . Figure [Tl] shows 
instead the reduced covariance matrix, Covn^didjd') = 
Cov/Aw{6)Aw{9'), obtained from the data with jack-knife 
(top panel) and analytically (bottom panel) They show a 
similar structure, with the jack-knife estimate being more 
noisy as expected. In addition at large angles {0 > 3°) there 
is a stronger covariance between separated angular bins in 
the jack-knife estimate probably due to systematics in the 
data. Nonetheless, as discussed in Sec. |4j this has no ma- 
jor impact in our study since the recovered best-fit models 
(and errors) derived using either jack-knife or analytical es- 
timates for the covariance are in broad agreement. In all, the 
underlying reason why the jack-knife and analytical error es- 
timates coincide is due to the fact that the model correlation 
functions are in good agreement with those measured in the 
data. This is the subject of the forthcoming sections. 



3.2 Star contamination 

Lastly, we study the star contamination of the sample as a 
function of the redshift. The quality of the star-galaxy sep- 
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Figure 11. Reduced covariance matrix estimated with 81 Jack- 
knife zones (top panel) or using the theoretical model in Crocce 
et al. (2011) (bottom panel). 



aration is a major concern in photometric catalogs since the 
broadband colors of stars can mimic those of galaxies and 
yield similar photometric redshift estimates, and important 
distortions in the clustering signal. 

We investigate the degree of contamination in our se- 
lected photometric sample performing the same selection, 
Eqs. ( 1|2|3 1, in the SDSS spectroscopic sample. In addition 
we only take the spectroscopic objects that overlap with our 
angular mask. For the bins where our analysis is performed 
we find f stars = 4 ± 1%. That is, a negligible dependence 
with redshift and a broad agreement with the residual con- 
tamination found in comparable clustering studies at these 
redshifts ( |Sawangwit et al.| p009| ; [Thomas et al.| ( |2011[ )). 

The next step is to estimate what is the impact of this 
contaminants in the large scale angular clustering signal 
since they introduce a density gradient through the galactic 
plane. Hence, we measure the correlation function of stars 
from the SDSS spectroscopic sample, relaxing the cut in 
magbO to have enough statistics. We also use the publicly 
available Tychos-2 star catalog ( |H0g et al.||2000a|b[ ) cut to 
the same selection and mask as our LRG sample to obtain 
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Figure 12. Angular correlation function of stars from the SDSS 
spectroscopic sample (error bars) and the Tychos-2 catalog (cir- 
cles) verifying our sample selection and mask. The correlation 
is well fit by Wstars,fit{9) = 0.0904 - 0.00313 9. Displayed error 
bars correspond to Poisson estimates, Aui = (1 -(- w) / \/Npairs 
(negligible for Tychos-2). 



a second estimate of the correlation of stars. Both determi- 
nations are in perfect agreement, as presented in Figure [121 
where the lines represent the correlation function for the 
SDSS sample, and the circles correspond to the Tychos-2 
catalog. This correlation is then included in the theoret- 
ical model for w{6) taking into account that LRGs and 
stars are uncorrelated populations, as (see also [Myers et al.[ 
(2006[)|Myers et al.[ ([2557|): 



Wobs,model{0 , z) = (1 — f stars) W gal .model {9 , z) + 

f star sW stars, fit{d) 



(8) 



where Wobs,modei is the model for the "observed" correla- 
tion function, Wgai^modei IS the model for the true correlation 
function of galaxies, Wstarsju is a simple fit to the measured 



correlation function for stars (see Figure 12 1. Notice that in 
Eq. (|8| we have removed the explicit dependence of the star 
fraction and correlation with redshift for simplicity. 



4 BIAS AND GROWTH OF STRUCTURE 

In this section we will employ the measured angular corre- 
lation function to place joint constraints in the growth rate 
of structure and bias of the LRG sample. 



4.1 Theoretical model 

We will use the theoretical model for the angular correlation 



function presented in Crocce et al. ( 2011 \. It was extensively 



tested against mock catalogs of photometric surveys with 
specifications similar to our present case. It was shown that 
the inclusion of redshift-space distortions and bias to linear 
order together with a model for non-linear matter clustering 
accurately reproduced the measured angular correlation in 
scales 8 > 0.5° for redshift bins centered at z ~ 0.5. We 



RSD and BAO with photometric LRGs from SDSS 9 



now recall some basic expressions of the model and refer the 



reader to Crocce et al. (20111 for further details 



The model angular correlation function is given by, 

Wgal,modetiO) ^ Jdzin{zi) j dz2n{z2) C^(l"l2) 



(9) 



where ri2 = ri2(zi, Z2, 9) is the separation of a pair of galax- 
ies at redshift zi and Z2 subtending an angle 6 with the 
observer. In turn, n{z) is the spectroscopic redshift distri- 
bution of the photometric sample under study. Notice that 
the inclusion of photo-z errors in the theory is solely through 
n{z) (e.i 



Budavari et al. (20031) 



Given the area and mean redshift of our sample we are 



allowed to make the plane-parallel approximation (see Rao 



canelli et al. (20101 and references therein for discussions of 
its validity). In this limit, the redshift space correlation is 



given by (Hamilton ( 1992 l) 



r(s,/.) =Co(s)Po(/x)+6(s)^'2(M)+e4(s)P4(M) (10) 

where Pi denote the standard Legendre polynomials, — 
r'^{zx)+r'^{z2) — 2r{z-i)r{z2) cos 9 and ^J. = {r{z2) — r{zi))/s, 
with r{z) being the comoving distance to redshift z. The £,e 
are the multi-poles of the spatial correlation 

Co(r) = (6'4-26//3 + /V5) K(r)] 

6(0 = (4&//3 + 4/77) [C(0-e'(0] 

C4(r) = (8/V35) [ar) + 5/2 ^'(r) - 7/2^(^)1 (11) 

with 5' = Sr-'^ i{x)x^dx and ^" = 5r~^' ^{x)x'^dx. 
Hence the angular correlation in Eq. ^ can be written as, 

w{9) = {b^ +2bf/3 + f/5)wo{9) + 

+ (46//3 + 4/V7) W2{9) + (8/V35) W4{e) (12) 

where we{9) are the bin projection of the functions in square 

urs 

bf and f, 

w{9) = b^woi9) + bf{2/3woi9)+4/3w2i9)) + 

+/71/5 wo(e) + A/7w2i9) + 8/35w4{9)). (13) 

The fact that each of these 3 terms have a different angu- 
lar dependence (through the different linear combinations 
of wt) makes it possible to constrain both b and / at the 



brackets in Eq. (111. Equation (12 1 can of course be re- 
arranged into three terms scaling as b^ 



same time. These terms are displayed in Fig. 13 (without 



the multiplicative factors involving b and /). The BAO is 
only present in the ^ = or "real-space" term determined 
by the bias. In real space the correlation becomes negative 
for scales larger than ~ 4° — 5° . The effect of RSD, in partic- 
ular of the cross term bf, is to make the correlation positive 
until larger scales, broadening the BAO feature (see also Fig. 
5 in [Crocce et al.| ( |2011[ )). Hence, at these scales the value 
of / can be degenerate with "excess-power" caused by sys- 
tematic effects that also make the correlation positive (e.g. 
extinction, star-galaxy separation or magnitude errors). 



4.2 Fits to growth and bias 

We will now use the measurements of the angular correlation 
function and the model described above to investigate the 
constraining power of the SDSS LRG photometric catalog 
onto the two-parameter space given by the velocity growth 
factor / and large-scale galaxy bias b. Recall however that 
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Figure 13. Angular dependence of the three additive terms mak- 
ing up the model angular correlation. They are proportional to 
b'^{z), b{z)f(z) and f'^{z) (in units of 0-8(2)) and are shown in 
solid red, dashed blue and dot-dashed black respectively. Their 
different shapes make it possible to constrain simultaneously b{z) 
and /(z). This figure corresponds to our photo-z bin [0.5 — 0.6]. 
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Figure 14. Two dimensional constrains in / X 0-8(2) and b{z) X 
0^(2). Dashed line corresponds to the 1 — a marginalized contour 
(Ax^ = 1) and solid line to the unmarginalized case (Ax^ = 2.3). 
In our fiducial cosmology 03(2 = 0.5) = 0.61, leading to best-fit 
values /(z = 0.5) = 0.87 and 6(2 = 0.5) = 1.84 



the multipoles ^ are arbitrarily normalized to, say, the am- 
plitude of fluctuations (js{z = 0) in spheres of 8ft~^ Mpc 
(e.g. Song & Percival (2009l). Thus, / and b are perfectly 
degenerate with this normalization and our two parameter 
space is in fact given by b{z)(jg,{z) and f{z)afi{z), where 
(jg,{z) = D{z)ofi{z = 0) and D{z) is the linear growth fac- 
tor. 

Notice that we will not make any assumption for the 
relationship between the velocity growth factor / and cos- 
mological parameters, but rather take / as a free-parameter. 
In particular we will not assume that the underlying theory 
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rcdshift bin 


b(z)as{z) 


f(z)as{z) 


dD/da 


0.15 - 0.30 


1.46 ±0.16 


0.49 ±0.10 


0.76 ±0.15 


0.30 - 0.40 


1.28 ±0.08 


0.42 ±0.06 


0.70 ± 0.10 


0.40 - 0.47 


1.46 ±0.16 


0.50 ±0.12 


0.90 ± 0.21 


0.50-0.60 


1.12 ±0.02 


0.53 ±0.42 


1.04 ± 0.81 



Table 1. Best-fit values for bias, velocity growth factor and 
growth rate. Top three rows correspond to the analysis of spec- 
troscopic SDSS LRG data as presented in Cabre & Gaztanaga 
(2009). Bottom row to our single redshift bin using the photo- 
metric LRG catalog. 



of gravity is General Relativity (where to a good approxi- 
mation / — Q^{z) and 7 — 0.55) because our ultimate goal 
is in fact understanding to what extent photometric data 
can help to constrain GR. 

This is different from other works in the literature where 
constraints from redshift space distortions for / (or /3 — f /b) 
are recast or combined with the ones for Qm assuming 



GR (e.g. Padmanabhan et al. (20071, Blake et al. (20071, 



Thomas et al. (20111, see also Sec. 4.4 1. Our approach cor- 
responds to the free growth model of [Samushia et ah] ( |2011[ ) . 
For the cosmological model we will assume a ACDM 



universe compatible with WMAP7 (Komatsu et al. (2011l) 
with = 0.272, Qde = 0.728, Qi, = 0.0456, = 0.963, 
h = 0.704 . In our approach these parameters determine the 
shape of the real-space matter correlation function and the 
distance-redshift relation. 

As shown in Sec. |2.3| splitting the data in multiple bins 
results in samples that are highly correlated. Thus, we de- 
cided to focus in a single redshift bin to present our main 
results and defer the study of narrow bins to Sec. |4.3[ in the 
context of robustness and consistency studies. 

Hence we concentrate on the bin [0.5 — 0.6]. This width 
is about twice the typical photometric error while the cen- 
ter of the bin makes this data uncorrelated with the SDSS 
LRG spectroscopic sample (that we refer to later on). The 
estimate for the true redshift distribution of galaxies in this 
bin is shown in Fig. [T] Using a spline-fit to it we computed 
the observed model correlation following Eqs. ( |8|9|12[ | sam- 
pling the two-parameter space given by fas — has ■ We then 
performed a standard minimization, where 

x\fas,bos) = A™(eOCov-/A«;(e,) (14) 

and Aw = ■Wobs,modei{f(7s,bas) — Wmeasured- To begin with 
we use the jack-knife covariance. The resulting best-fit values 
and 1 — a errors are listed in Table [l] (to convert to / and b 
one can use that (78(0.55) — 0.611 for our fiducial cosmology 
iiaaiz = 0) = 0.8). The fit yields xli^n = 26.67 for 28 bins in 
9 in the range [0.6° — 6°] and two fitting parameters. Hence 
the quality of the fit is very good. The joint error distribution 
is displayed in Fig. [14] and shows that these two parameters 
are not degenerate. The significance of RSD in our data is 
~ 1.26 a" (i.e. the degree by which the recovered value for / 
is away from zero). 

Alternatively one can assume that the velocity growth 
rate is given by General Relativity, i.e. f{z)as{z) — 



(J8(2:)dlnD(2:)/dlna = 0.45 (for our fiducial cosmology and 
(Jsiz = 0) = 0.8), and fit only for the overall amplitude of 
clustering. In this case we recover a similar Xrnin{~ 26.7) 
and has ~ 1.125±0.017. Hence the fit does not vary appre- 
ciable, showing that agreement with standard ACDM is also 
very good and that the sensitivity of our data to redshift- 
distortions is weak, as expected. 

The results so far were obtained assuming 4% star 
contamination (as discussed in Sec. 3.2 and Eq. (|8|). If 
we now assume that our sample is perfectly clean of stars 
if stars = 0) we get bas = 1.08 ±0.02 and fas = 0.66 ±0.39 
and a Xmin ~ 26.87. This model is hardly differentiable from 
the case with f stars ~ 4% (the x^ becomes only negligible 
worse), see Fig. 15 Hence stellar contamination can mimick 
the effect of RSD and become a worrisome source of sys- 
tematic effect in future and better data. In the present case 
the agreement between the best-fit value of / and the "GR" 
value degrades only slightly when fstar = 0, with differences 
well within the errors. We conclude that the star contami- 
nation in our sample is sufficiently under control and does 
not drive our results. 

Figure [15] shows the resulting best-fit models discussed 
above against the measured angular correlation function. 
In all cases the model matches the data quite well in all 
the angular range, with no signal of excess-clustering on 
large scales. This is to some extent at variance with the 



recent study by Thomas et al. (20111 (see also the follow 
up Thomas et al. (2010l) who finds the amplitude of the 
angular power spectra of photometric LRGs to have a 2- 
a excess clustering away from the ACDM prediction at the 
lowest multipoles. This result was obtained using the MegaZ 
catalog over three redshift bins in the range [0.4-0.65]. The 
MegaZ catalog is also build upon the DR7 SDSS photometric 
catalog. However the LRG selection criteria in MegaZ is dif- 
ferent from the one in this paper, most notably by the magni- 
tude cuts (idcv < 19.8 in MegaZ compared to petror < 21 in 
our case) but also in the color and mag^o cuts intended to 
isolate stars from galaxies. The photo-z code used for MegaZ 
is also different from ours, although both arc based in the 
ANNz code of |Collister fc Lahav] ( |2004[ ). [Sawangwit et al. 
( 2009 ) also finds an excess clustering when searching for the 
baryon acoustic scale, but this result is difficult to compare 
with ours as they use an over simplified method to stack 
the signal from wider redshift bins calibrated with different 
spectroscopic samples. 

The constraints recovered so far, albeit large, match the 



study of Ross et al. (20111 that forecast how well redshift 
distortions can potentially be measured with galaxies se- 
lected from a photometric survey. This match can be clearly 
seen from their Fig. 8 after scaling the results for a sample 
of unbiased galaxies binned around z = 0.5 by our bias fac- 
tor & ~ 2 (since A{fas) oc b). For our redshift bin width 
the forecast yields A (/erg) ~ 0.4 — 0.5, in very good agree- 
ment with the results we obtain with the actual data (e.g. 
see Table [T|. This implies that systematics, photo-z errors, 



The MegaZ selection is done to match the magnitude cut in 
2SLAQ from where the photo-z calibration is obtained. We find 
that the angular correlation of galaxies in MegaZ matches the 
clustering amplitude in our sample if we cut in petror < 20.7. 
Hence MegaZ is slightly brighter than our sample. 
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Figure 15. Angular correlation function in our central bin of width 0.1 centered at z = 0.55. Solid red line is our best-fit model including 
the effect of star contamination {/star = 4%). Solid blue line is the corresponding best-fit model if fstar = 0. The values for f and b are 
given in Table ^ For reference we include with a dashed black line a WMAP7 ACDM model assuming General Relativity, tha t is with / 
set to n^^^ {fstai- = in this case). The inset panel zooms in the region where the baryon acoustic peak is located (see Fig. |l8| for the 
BAO significance). Notably the different models match the data very well in all the range of scales. 



selection and modeling can be controlled sufficiently well in 
present photometric data to yield expected constraints of 
redshift distortions. Near future data, provided with a bet- 
ter handle on systematics due to stars and photo-z, should 
be able to complement shallower spectroscopic surveys in 
placing bounds to the growth of structure in the universe. 



4.3 Robustness and impact of systematics effects 

In this section we investigate different components of our 
analysis that could potentially change our results. 

Analytical Covariance Matrix: We now investigate how 
our results change when we use the analytical error estimate 
discussed in Sec. ( |3.1[ ). This estimate accounts for statistical 
and shot-noise errors but does not account for systematic er- 
rors introduced by, say, stars (at least in the way presented 
in |Crocce et al.| ( |2011 1). Hence we compute the Ce spectra, 
and CovT/i, in Eq. ([7i using the best-fit values for / and b 
corresponding to the model with f stars = discussed above. 
Using this covariance the best-fit values to the angular cor- 
relation in the bin [0.5 — 0.6] change to bas ~ 1.13 ± 0.02 
and fas = 0.35 ± 0.54 while the Xmin degrades to 32 (for a 
model with /stars ~ 4%). Hence the change in the recovered 
best-fit values is within half a compared to the results shown 



in Table [T] In turn, the error in has is unchanged while that 
in fas increases by about 25%. 

The resulting model matches the data as well as the 
one derived using jack-knife errors and the goodness of the 
fit are comparable. The only difference is in the resulting 
error on the velocity growth factor. This may be due to 
the structure of the covariance on large separations that 
may be affected by some systematics not captured by the 
theoretical estimate. 

Redshift Distribution: One important but difficult to es- 
timate component in clustering analysis of photometric data 
is the distribution of galaxies in true redshift. In Sec. |2.3| 
we studied the degree of uncertainty left in the estimate of 
N{z), about 1% in the peak position and 9% in the width. 
To investigate the impact of this systematic in our analysis 
we computed the model correlation and found best-fit 
parameters assuming N{z) is a Gaussian distribution with 
H — 0.549 and a^ = 0.062 (model-1), that are the best-fit 
values to N{z) in Fig. [7] We then recomputed the model by 
increasing n by 1% and decreasing a^ by 9% (model-2), and 
subsequently did the best-fit analysis. The of these two 
models is almost equal, changing by ~ 1%. The resulting 
best-fit values for bas decrease from 1.12 ± 0.02 (for 
model-1) to 1.08 ±0.02 (for model-2) while from 0.76 ±0.57 
to 0.7 ± 0.5 for fas- This is because a narrower N{z) 
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Figure 16. Angular correlation function in "narrow" bins of 
width comparable to the underlying photometric error, as indi- 
cated in the panels. Solid red lines are the best-fit models taking 
into account the ~ 4% star contamination in our sample. Dashed 
line correspond to a standard ACDM model with fstar = and 
growth given by General Relativity, i.e. f{z) = Qrniz)^'^^ ■ 



implies less bin proyection what leads to a higher amphtude 
angular correlation and slightly more sensitivity to redshift 



distortions ( |Nock et a.l\ ( |2010[ ) ; [Crocce et aL] ( |2011| ) . These 
differences are only marginal given the overall large er- 
rors provided by present photometric data. However future 
surveys will probably need more accurate estimates of N{z). 

Narrower redshift bins: As mentioned earlier we have 
also considered splitting the data into 3 narrower bins of 
width similar to the typical photo-z error : [0.45 — 0.5], 
[0.50 - 0.55] and [0.55 - 0.60]. The galaxy redshift distribu- 
tions for these bins are shown in Fig. |8] while the measured 
correlations and best-fit models are displayed in Fig. |16[ 

The best-fit values for the bias in each bin decreases 



slightly with increasing redshift: bag = 1.26, 1.21 and 1.1 
respectively (with 2% error). 

In turn the best-fit values for /erg are 1.14±0.57, 0.024± 
0.53 and 1.39 ± 0.46 respectively (assuming fstar = 4%). 
Hence we see some spread in the recovered values for the 
velocity growh rate. Compared to the corresponding values 
in GR they still agree at ~ 1 — 2(t. 

As discussed before a bin width smaller or comparable 
to the photo-z is subject to large bin to bin migration (in 
other words, the estimate of N{z) itself is more sensitive 
to photo-z unkowns). Therefore we expect to recover more 
robust results in bins larger than the underlying photo-z. 



4.4 Implications for the growth rate and 
comparison with spectroscopic studies 

We now put our results in context of similar ones derived 
using spectroscopic data. In Table [l] we show constraints in 
b{z)as{z) and f{z)D{z) as obtained by |Cabre &: Gaztanaga| 
( 20091) using spectroscopic LRGs from SDSS. Table 2 in 



Cabre & Gaztanaga ( 2009 1 list the best fit values of the 
clustering amplitude Amp = b{z)as and Qm- This Table 
refers to (Tg at z = because the best fit value of was 
used to estimate the linear growth D{z) at the correspond- 
ing redshift according to standard cosmological equations in 
General Relativity (GR). Here we do not want to assume GR 
or any other relation between Om and D{z) and we therefore 
scale these amplitudes back to the original data amplitudes, 
i.e. crs{z), by multiplying Amp by the best fit value of D{z) 
in GR. The resulting amplitudes are listed here as b{z)as{z) 
in Table 1. We can then find an estimate of f{z)crs{z) by 
just multiplying these b{z)ag{z) estimates with the values 
of /3 listed in Table 1 of Cabre & Gaztanaga ( 2009 1 for the 
corresponding samples. 

One can turn these values, and those a,t z = 0.55 from 
the photometric data, into estimates for the linear growth 
rate of structure as follows, 

dD D{z) ^ + ^ ft , ( ^ 

where we assume our fiducial value crg(O) = 0.8 (consistent 
with, 



e.g.. 



Tinker et al. (2011l). Results are listed in Ta- 



ble [T] and displayed altogether in Fig. |17| Provided with 
them one give constraints in the growth index assuming the 
7-parameterization of the growth, for which / = Q.m{zy 
and 7 is a scale and redshift independent constant, |Linder| 
(20051). In this case, 



dP 
da 



(l + z)D(z)Qm{zy 



and D{z) — exp (— J^^ Q.rn{x)'^ / (1 + x)dx). This yields 
7 = 0.54 ±0.17 



(16) 



(17) 



for the combination of spectroscopic and photometric data 
given in the fourth column of Table [l] 

Before moving on we note that novel constraints on the 
growth rate up to 2 = 0.9 were very recently reported by 
the WiggleZ survey using blue galaxies instead of LRGs 



(Blake et al. (2011 1). At 2 < 0.5 they improve to some ex- 
tent over the ones we used, e.g. A(/crg) = 0.07 at z = 0.2 
and A(/a8) = 0.04 at 2 = 0.4 (see Tablejl]). At 2 > 0.5 they 
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Figure 17. Linear growth rate of structure from LRG spectro- 
scopic data in the range [0.15-0.47], as presented in Cabre & 
Gaztariaga (2009), and from our analysis of photometric data 
at 2 = 0.55 (assuming a-g(O) = 0.8). These data leads to 
7 = 0.54 it 0.17 in a model where / = Qrn{z)'^ ■ 



are considerably tighter than those we derive with imagining 
data, as expected. Yet, for concreteness we have decided to 
present our results in terms of SDSS LRG clustering (either 
spectroscopic or photometric) as the focus of the paper is to 
demonstrate that imaging data is able to yield constraints 
in RSD. 



5 BARYON ACOUSTIC OSCILLATIONS 

One of the most important probes of the accelerated cos- 
mic expansion is the existence of an excess clustering im- 
printed at a well defined comoving length scale of about 
~ 100/i~^Mpc in the correlation of galaxies and ~ 1° in 
the one of CMB photons. It originates in the coupling of the 
baryon-photon plasma prior to recombination and hence can 
be very well determined with CMB data. Provided with this 
estimate it can then be measured with galaxy data and used 
to constrain the distance-redshift relation in the local uni- 
verse, that in turn is sensitive to the nature of dark-energy 
(and/or the appropriate law of gravity). The main obstacle 
of this pathway is that the excess clustering signal in galax- 
ies represents only ~ 1% over that of a random distribution. 
Nonetheless several surveys are dedicated or include BAO 
in their scientific plans. 

The BAO signature have been extensively studied in 
the clustering of spectroscopic data, in particular using 
LRGs since they span the largest volume possible compared 



to other galaxy types (e.g. Eisenstein et al. 
(2006);'Percival ot al.'(2007')''Gaztanaga ct al 



2005h; Hiitsi 



|& Gaztaiiaga (2009); Gaztariaga ct al. 



2009 j); [Cabre 
2009): Kazin et al. 



(20101; Percival et al.| (2010|). The advantage in this case is 
that one has three dimensional information to sample this 
1% excess of pairs. The disadvantage is that the sampling of 
the over-density field is much more poorer with spectra and 



is eventually limited to lower redshifts when compared to 
photometric catalogs (which on the converse only yield pro- 
jected quantities with lower significance). Hence it is very 
important to investigate what evidence for this signature is 
already present in our data given the number of such photo- 
metric surveys already undergoing or planned for the near- 
future. 

The use of photometric data to investigate BAO has 
been relegated to some extent, probably due to the lower 
signal-to-noise, the impact of systematics and the quality 



of the photo-z. In Blake et al. (20071 and Padmanabhan 



et al. ( 2007 1 this signature was studied stacking photomet- 
ric redshift bins in order to reconstruct the 3D spectrum. 
Both studies find evidence for BAO at < 3ct. In turn, only 



Sawangwit et al. ( 2009 1 (to our knowledge) have explored 
the possibility of locating the signature in configuration 
space, with ambiguous results (e.g. too high in amplitude 
compared to ACDM expectations). 

We will now investigate how well our measurements 
agree with the shape of the model correlation including 
or excluding the effect of baryons. In a separate paper we 
present a detailed analysis of the significance and implica- 
tion of the baryon acoustic peak imprinted in the angular 



clustering of our sample, see Carnero et al. (20111 



Figure [18] shows the measured angular correlation in the 
photo-z bin [0.5 — 0.6] in a way that highlights the shape 
of w{6) at large angular scales. We also display the stan- 



dard Eisenstein & Hu ( 1998 1 model prediction including the 
baryons effect (solid blue) and excluding it, i.e. with the wig- 
gles smoothed out (solid red). Figure 18 shows a quite clear 
bump in the correlation function at 6* ~ 4°, in very good 
agreement with the model in both shape and location. For 
this figure we have not attempted to fit for redshift distor- 



tion parameter f{z) as in Sec. 4.2 but rather have assumed 
GR {f{z = 0.55) = 0.74) and fit for the overaU bias only 
(with fixed cosmological parameters as given in Sec. 4.21. 
We note that the "BAO" model correlation obtained with 



the approximation of Eisenstein & Hu ( 1998 1 is very similar 



the one derived with CAMB in Sec. 14.21 

To estimate the statistical significance of this feature 
we compare the two best-fit models, with and without BAO. 
The difference in their yields 5.4 suggesting that the BAO 
model is preferred with a significance ~ 2.3 a (~ 98% con- 
fidence). Notice that both models have the same number 
of degrees of freedom since only the bias is fit. This re- 
sult is un-changed if we allow for 4% star contamination 
in both models and is in perfect agreement with the find- 



ing of I Carnero et al. (2011 \ using a completely independent 
analysis method. 

Although the significance is low the good agreement 
with ACDM expectations is very encouraging for future pho- 
tometric campaign that will achieve better error bars thanks 
to improved photo-z and survey depth. In addition this is, to 
our knowledge, the first time the BAO bump is clearly de- 
picted in agreement with ACDM predictions in the angular 
correlation of photometric LRG data at this high redshift. 



6 CONCLUSIONS 

In this paper we used the angular correlation function of 
the imaging sample of luminous red galaxies in the DR7 
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be taken as a validation of the forecast for a SDSS like-case 
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as presented in Ross et al. (20111 and a proof-of-concept 
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Figure 18. Baryon Acoustic Feature at z = 0.55. The plot shows 
the Eisenstein & Hu (1998) model with and without baryons. The 
data and model agree quite well in both amplitude and shape of 
the acoustic bump. The BAO feature is detected here with a 
98% confidence level (2.3 cr significance). Dashed line corresponds 
to the no-wiggle model when RSD are ignored. This shows the 
importance of both BAO and RSD to account for observations. 



of the SDSS as a testing ground for measuring the growth 
rate of expansion through redshift space distortions using 
photometric data. In addition we investigated the evidence 
for the baryon acoustic feature in the angular correlation. 

We put a strong emphasis in the selection of the galaxy 
sample and overall robustness against several systematic ef- 
fects, such as magnitude errors, bad extinction zones, photo- 
z outliers and more. We paid particular attention to the 
impact of stellar contamination in the angular clustering 
measurements. On the one hand we minimized such con- 
tamination introducing a cut in surface brightness, magso- 
On the other we estimated the distortion that the residual 
contamination introduces in the correlation measurements. 

Our measured correlation, in the range [0.50 — 0.60] is 
in good agreement with expectations from standard ACDM. 
In particular it shows no excess clustering on the largest 
scales, contrary to other works in the literature. This is a en- 
couraging proof that systematic effects in photometric data 
can be controlled sufficiently well to use them as a cosmo- 
logical tool ([Blake et aL] ( |2007[ ); [Padmanabhan et al.| ( [20071 ); 
Thomas et al. (2011 1). The distortions introduced by the in- 
trinsic clustering of contaminating stars does not change our 
results but might become a worrisome source of systematic 
biases in future surveys with smaller error bars. A similar 
conclusion was reached in regards to the estimation of the 
true redshift distribution of the sample. 

Indeed, we found that redshift space distortions can be 
measured using photometric data, albeit with large error 
bars due to the high-bias of the sample and the poor photo- 
z error. Our results are in very good agreement with the 



recent forecast by Ross et al. (20111. Hence, this paper can 



of the promising expectations for upcoming photo-z surveys 
such as DES, Euclid and PannStars. 

In addition we found quite a strong evidence for the 
baryon acoustic peak in the measured angular correlation, 
something not observed before. The shape, amplitude and 
location of the BAO feature is in very good agreement with 
ACDM expectations yielding a ~ 2.3 ct significance over a 
model without BAO. In a separate work, [Carnero et al.| 
( |201l| , we discuss this detection of BAO and its cosmo- 
logical implication using an independent analysis from the 
one presented here. 

In all, our results strengthen the expectations on the 
ability of future photometric surveys to compete and/or 
complement spectroscopic data, as well as to serve to other 
approaches such us weak lensing or supernovae, in the quest 
to understand the nature of cosmic acceleration. 
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APPENDIX A: EXCESS POWER FOR Z > 0.6 

One frequent issue when analyzing the clustering of LRGs is 
the fact that at the largest scales (e.g. BAO) the amplitude 
of clustering appears generally high when compared to stan- 



dard ACDM models (Eisenstein et al. (20051; Blake et al. 


(20071; Padmanabhan et al. 




20071; Thomas et al. (2010 


2011 1; Sawangwit et al. ( 2009 




Okumura et al. ( 2008 1; Kazin 


et al. 


( 


20101; Cabre & Gaztafiaga (20091; Sanchez et al. 


(2009 


1; 


Samushia et al. (2011 1). The significance of this mis- 



match is however uncertain since these scales are expected 
to be the most sensitive ones to different systematic uncer- 
tainties as well as the ones with largest statistical variance. 

Within the context of photometric data as in our work, 
Thomas et al. ( 2011 2010 1 recently found ~ 2a excess in the 
lowest multipoles of the angular power spectrum of LRGs in 
the MegaZ catalog of DR7 SDSS. As discussed in Sec.[4]we 
do not find such a discrepancy in our catalog. More so if one 
recalls that in configuration space data points are expected 
to co-variate to some level. Potentially much more worrisome 
is their finding of ~ 4it excess for the bin [0.6 — 0.65] (see 



also Blake et al. (20071). Thomas et al. (20101 performed a 



series of checks for systematic errors but none was conclu- 
sively the source of such an effect what led them to speculate 
with the possibility that this could be due to the imprint of 
"new physics" , such as primordial non-Gaussianities, modi- 
fications of gravity or clustering dark energy. 

In Fourier space this excess may not be prejudicial be- 
cause it only affect few low-^ multipoles that can be cut-out 
of the analysis (they have the lowest signal-to- noise anyway). 
When translated to configuration-space this impacts a broad 
range of scales. Hence it is interesting to see how does the 
clustering in configuration space looks like ai z ^ 0.6 to 



complement the study of Thomas et al. (2011 20101 



In Fig. |A1| we show the measured angular correlation 
function for our photometric bin [0.6. — 0.65]. Evidently the 
clustering signal does not only show an excess power at BAO 
scales but is in fact anomalous at all scales (except perhaps 
6 < 0.5°). On the one hand the number of objects in this 
bin is small (~ 30, 000) and dominates the error budget 
(shot-noise). On the other hand the correlation is very sen- 
sitive to various systematic uncertainties. This is reflected in 
Fig. |Al| : the solid magenta line corresponds to the measured 
correlation when no cut in galactic latitude of the sample or 
photo-z quality is imposed. Blue solid line is the result when 
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Figure Al. Excess power beyond z = 0.6 and sensitivity to sys- 
tematic effects. Top magenta line corresponds to a mask with no 
cut in galactic latitude or photo-z quality. Middle blue line cor- 
respond to a cut in photo-z quality as discussed in Sec. |2.3| Red 
solid line corresponds, in addition, to a mask excluding galac- 
tic latitude below 25°, see Sec. |2.l| For reference we include the 
measured correlation in the bin [0.5-0.6]. The correlation seems 
anomalous showing an excess power on all scales and a large sen- 
sitivity to various systematic effects. 



objects with "bad" photo-z are discarded from the sample 
(as discussed in Sec. 2.31. Lastly, solid red line is the result 
when the mask is reduced by leaving out low galactic lat- 
itudes b < 25° (to avoid star contamination, see Sec. |2.1[ ). 
One could continue with more stringent constrains but the 
signal does not approach the one at lowest bins (shown with 
black symbols). Displayed error bars in this figure were ob- 
tained with jack-knife estimate. 

This result signals to unknown systematic uncertainties 
as the most probable cause for the anomalous shape. One 
possibility could be an incomplete treatment of the distor- 
tions introduced by star contamination. This is generally 
focused in reducing the fraction of stars to the minimum 
possible but not to account of the residual clustering sig- 
nal, as in Sec. |3.2| The contamination by stars results in a 
change of the mean density across the survey area (and to- 
wards the galactic plane) . This would impact only the low-^ 
spectra that encodes the mean density information but a 
broad range of scales in the angular correlation (as shown 



in Fig. 12 1. However testing this in Fourier space is beyond 



the scope of this paper and will be discussed elsewhere. 



